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30% of the DNA in E. coli bacteria is covered by proteins. Such high degree of crowding affect 
the dynamics of generic biological processes (e.g. gene regulation, DNA repair, protein diffusion 
etc.) in ways that are not yet fully understood. In this paper, we theoretically address the diffusion 
constant of a tracer particle in a one dimensional system surrounded by impenetrable crowder 
particles. While the tracer particle always stays on the lattice, crowder particles may unbind to a 
surrounding bulk and rebind at another or the same location. In this scenario we determine how the 
long time diffusion constant V (after many unbinding events) depends on (i) the unbinding rate of 
crowder particles k 0 g, and (ii) crowder particle line density p, from simulations (Gillespie algorithm) 
and analytical calculations. For small k a g, we find V ~ k 0 g/p 2 when crowder particles are immobile 
on the line, and V ~ \JDk 0 a/p when they are diffusing; D is the free particle diffusion constant. For 
large fc 0 fi, we find agreement with mean-field results which do not depend on k 0 g. From literature 
values of k a g and D, we show that the small fc 0 ff-limit is relevant for in vivo protein diffusion on a 
crowded DNA. Our results applies to single-molecule tracking experiments. 


I. INTRODUCTION 

Few doubt that molecular crowding has severe conse¬ 
quences for dynamical processes [T] . Interesting examples 
are living cells where macromolecular concentrations are 
large. Take the E. coli bacterium as an example. There, 
the concentration of proteins and RNA is about 300-400 
mg/ml 0 which is 30-40 times larger than common test 
tube conditions |3j. There is overwhelming evidence that 
this level of crowding influences important biological pro¬ 
cesses such as gene regulation [3], enzymatic activity [5], 
protein folding 0 0 , and diffusion of macromolecules 
BM- In order to get a complete picture of the in vivo 
dynamics we must increase our understanding of the role 
of crowding, and recent experimental developments pro¬ 
vide the means to do it. 

In recent years, researchers have beat the diffraction 
limit and turned optical microscopy into ’nanoscopy’. 
Today’s microscopy methods (e.g. STED, STORM and 
FIONA) [TO] does not only allow us to image nanometer¬ 
sized biological structures, but recent improvements 
mm also permit tracking fluorescently labelled pro¬ 
teins at the biologically relevant millisecond-scale. This 
is anticipated to shed new light on biological processes, as 
well as increase our understanding of particle transport 
in engineered nano-fluidic systems mm- In order to 
properly interpret those type of experiments in vivo, we 
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need new theoretical and computational models in terms 
of physical properties of the intracellular space, the cy¬ 
toplasm. 

The cytoplasm is cramped with macromolecules and 
we are interested how this influences diffusion-controlled 
processes, a key component in many cellular functions 
(e.g. gene regulation). While our results are new, as¬ 
pects of this problem has been studied theoretically be¬ 
fore. For example, mm investigate diffusion in the 
three dimensional cytoplasm and in gels, whereas [T01 - I22| 
focus on the sub-diffusive motion seen in single-molecule 
experiments. Crowding is also important for DNA search 
processes where a searcher combines one and three di¬ 
mensional diffusion to quickly find its target, so called 
facilitated diffusion. Facilitated diffusion under crowd¬ 
ing is addressed in 0122 ] which resemble this paper but 
we ask different questions: we calculate the diffusion con¬ 
stant of a tracer particle in terms of key properties of sur¬ 
rounding crowder particles rather than focusing on mean 
target finding times. 

Much inspiration to this work comes from DNA bind¬ 
ing proteins. Of particular interest is repair proteins 
(MutS and homologs) whose residence time on the DNA 
can be very long (~ 10 min [23j ). We are also inspired by 
transcription factors, the family of gene regulatory pro¬ 
teins. The yeast regulatory proteins LexA and Gal4 can 
stay bound to their regulatory sites for several minutes 
in vitro (LexA ~ 5 min and Gal4 > 30 min) [25] , but sur¬ 
prisingly, this number can be reduced up to 1000 times in 
vivo. Both classes of proteins have the ability to diffuse 
along the DNA, unbind to the three dimensional intra- 
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FIG. 1: Schematic illustration of our model. All particles are 
diffusing with rate ko on a one dimensional lattice with lattice 
spacing a. The crowder particles (green) may also unbind and 
rebind to a random, or the same, lattice site with rates k 0 s 
and fc on , respectively. The tagged particle (orange) cannot 
leave the line (fc on = fcoff = 0). 


cellular space, diffuse in space, and rebind to the DNA. 
We are interested in how the dynamics of those proteins 
change under crowding. 

In order to better understand the role of crowding, 
we introduce a theoretical model where particles diffuse 
on a one dimensional lattice where two particles cannot 
occupy the same site (Fig. [l] ). They diffuse with rate 
kp (same for all particles) and may unbind and rebind 
to the lattice with rates k Q g and fc on , respectively. These 
rates are tuned so that the average particle line density 
is constant at 10-20% which is not too far from in vivo 
conditions (30% of the DNA in E. coli is covered by 
proteins). The unbinding rate for the tracer particle 
is set to zero similar to the long-lived protein-DNA 
complexes described above. Now we ask: 

What is the long time diffusion constant of a tracer 
particle in such a crowded quasi one dimensional system? 


II. THE MODEL 

Our model has been used and explained elsewhere [34], 
but for completeness we summarise it briefly below. Con¬ 
sider a one dimensional lattice on which crowder particles 
(assumed identical) and the tracer particle diffuse (Fig. 
[Tj). The crowder particles can diffuse, unbind and rebind 
to the lattice. Rebinding occurs in two ways. Either to 
a random unoccupied lattice site (chosen with uniform 
probability), or to the exact same location. Both rebind¬ 
ing modes has been used to model transcription factor 
dynamics on DNA [.‘351 [5S] , and we will therefore con¬ 
sider both. The lattice constant is denoted a, and the 
diffusion rate kg, is assumed equal in both directions and 
for all particles. Double occupancy is forbidden and a 
particle cannot overtake a flanking neighbor (single-file 
condition). Binding and unbinding dynamics of Crow¬ 
ders are characterised by the rates k on and fc 0 ff, which 
are chosen such that the particle line density is in equi¬ 
librium with the bulk, thereby keeping the average filling 
fraction is constant. In our simulations we keep it at 
10-20%. We implemented the model using the Gillespie 
algorithm. See Appendix [C| for details. 

III. ANALYTICAL ESTIMATES FOR THE 
LONG TIME DIFFUSION CONSTANT 

Here we provide analytical estimates to corroborate 
and better understand the numerical results in the next 
section. We are mainly interested in the long time diffu¬ 
sion V constant for the tracer particle, defined as 

(x 2 (t)) ~ 2 Vt (1) 


We answer this question numerically using stochastic 
simulations (Gillespie algorithm), corroborated with an¬ 
alytical results. The main results are Figs. [4]-|6] where we 
show how the diffusion constant changes as a function of 
our main parameter k D g. Those results are applicable to 
single molecule tracking experiments }26j . 

The unbinding rate k 0 g interpolates between two well 
studied limits, (i) When fc off is large (compared to kff), 
the tracer’s mobility is only weakly lowered and diffuses 
close to as if it was free, (ii) When k 0 g —>• 0, the particles 
diffuse with unchanged order in a single file. Single-file 
diffusion is well studied [27H53] where the most famous 
result is that the mean squared displacement of a tracer 
tracer particle is proportional to yft rather t (t is time) 
which signatures non-markovian dynamics. 

This paper is organised as follows. In Sec. II, we out¬ 
line briefly the details of our model. Before showing the 
results in Sec. IV, we provide analytical estimates of the 
diffusion constant in Sec. Ill, based on a theoretical cal¬ 
culation found in Appendix [A] In Sec. Ill we also briefly 
review the dynamics of the model at short, intermediate 
and long times. We close by a few concluding remarks in 
Sec. V. 


where ( x 2 (t )) is the ensemble averaged mean squared dis¬ 
placement (MSD), and t is time. Notably, V is in general 
not equal to the bare, or free particle, diffusion constant 

D = a 2 k D . (2) 


It is a non-trivial function of k 0 g and p. To better under¬ 
stand what we mean by long time, we describe in sub¬ 
section [nllD the dynamics leading up to Eq. (JT|). But 
first we summarise our main analytical findings from Ap¬ 
pendix [A] which we in Sec. |IV| compare to simulations. 


A. Long time, small k 0 g behaviour 

To simplify matters, we start by assuming that the 
crowder particles sit equidistantly on the line with den¬ 
sity p, unable to diffuse (ko = 0), and rebind to the site 
from which they unbound. In this situation, the tracer 
moves back and fourth between its flanking neighbours 
and can only move past them if one of them unbinds. 
The average time until this happens is proportional to 
In point of view of the tracer this process is a 
random walk on an effective, or coarse grained, lattice 
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with spacing and jump rate proportional to 1/p and fc 0 s, 
respectively. From this we expect that V ~ k Q s /p 2 , and 
a more elaborate calculation shows that 


3(1 - ap)k 0 g 
4p 2 


( 3 ) 


When crowder particles rebind to a random location 
rather than to the same site, the distance between two 
neighbouring particles fluctuate even though the average 
density is fixed. This leads to a larger effective lattice 
spacing, and a larger V compared to Eq. ®>: 

(2-ap)(i-ap) 

V ~ - — 2 - k oS . (4) 

When crowder particles also diffuse (kn ^ 0) the dis¬ 
tance between nearest-neighbours becomes difficult to 
define. We estimate the coarse grained lattice constant 
as the length the tracer particle explores during a time 
intervall proportional to 1 /k 0 s- This leads to 


V * (1 ~ J™*, (5) 

p V 27r 

which has different fc 0 ff-scaling than before. Equations 
© © constitute our main analytical results. 


B. Long time, large k 0 n behaviour 

When kjj > fc 0 g, crowder particles frequently unbind 
and rebind to the lattice and the no-passing condition 
is effectively violated. But, crowder particles still hinder 
the tracer thereby decreasing the diffusion rate. Imagine 
that the jump rate for a single particle to a neighbouring 
site on an otherwise empty lattice is kjj, or D = a 2 ko- 
Then, when crowder particles are around, some of the 
jumps are canceled because the target lattice site may be 
occupied. In that situation the jump rate is reduced by 
the probability that the target lattice is unoccupied. For 
very large k 0 g this probability is simply 1 — ap , therefore 

V ~ D{1 - pa) (6) 

This mean field result has been obtained before nama 
[55] . and as k Q s/ku is close to or smaller than unity, cor¬ 
rections to this formula becomes prominent (see Fig. [5]). 


C. Interpolation formula for V 


Based on the expressions above, we propose a simple 
formula for V valid for all k 0 g: 


1 

V 


+ 


T^small fc Q ff -^large feoff 


( 7 ) 


Here fc^iarge fc 0 ff Eq. © whereas 2? sma n fc off is one of 
Eqs. j3])-|5]) depending on the case under study. Equa¬ 
tion (0) is appealingly simple and captures properly the 


small and large k 0 g limits, but it should not be viewed 
as more than a candidate expression for T>. We have not 
made a systematic attempt to find the best form of U, 
and leave it for future research. 


D. How is the long time asymptotics [Eq. ([l])] 
approached? 


Here we clarify the meaning of short, intermediate and 
long times within of our model. To keep the discussion 
simple, we consider fc Q ff, k on , kjj, and p as constant. See 
also Fig. 0 which shows (x 2 (t)) as a function of time, 
where all relevant regimes are present. 

At most we have three regimes of different behaviour. 
These are separated by the average residence time of the 
crowder particles T 0 ff, and the average collision time r co u, 
which is the time it takes for a particle to diffuse across 
the average nearest neighbour distance 1/p: 


Tcoll 


1 _ 1 

P 2 D ’ T ° ff fc off 


( 8 ) 


Let us assume that there is a clear separation between 
these timescales and that r co ii <C r 0 ff and that kp is the 
fastest rate in the system, 1/kjj <C r co u. In the first 
regime, t <C t co ii, the tracer diffuse as if it was free, since 
it has not yet collided with its nearest neighbours. This 
means that the tracer’s MSD is {x 2 (t)) = 2 Dt. In the 
second regime, r co n < i < r 0 ff, many particle collisions 
have taken place but particles diffuse with maintained 
order since they are unable to pass each other. This is 
the single-file diffusion regime which is characterised by 
Harris’ law ( x 2 (t )) oc \jDt/p 2 [27]. Here, memory ef¬ 
fects dominate and is the very reason to the sub diffusive 
behaviour. In the third regime, t 3> T 0 ff, particles start 
unbinding from the lattice which effectively violates the 
no-passing condition. In this regime we expect diffusive 
behaviour again (x 2 (t)) ~ t, but with a diffusion con¬ 
stant different from D , denoted by V [see Eq. ©]• This 
is the one we wish to calculate, in particular in terms of 
our key parameter fc Q ff. Note that the second regime can 
be erased completely if we lower r Q ff such that r 0 ff « r co ii 
(or smaller). Similarly, the third regime is absent if un¬ 
binding is not allowed, i.e. fc Q ff = 0 (or r 0 g = oo). In 
most of our simulations, diffusion is the fastest process 
in the system which also is the likely scenario a biological 
cell (see Sec. [0). To sum up, 

r 2 Dt, t < T CO ii 

(x 2 {t)) ~ < (1 - pa)/p x y/4/DtJ/r, Tcoll < A f < A T 0 ff 

[ 2 T>t, t » Toff 

(9) 


IV. RESULTS 

In this section we present results from stochastic sim¬ 
ulations of the model outlined in Sec. |TIJ together with 
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of our theoretical findings from Sec. |III| The simulation 
details can be found in Appendix [C] First, we show the 
tracer particle’s MSD as a function of time, from which 
we extract the long time diffusion constant V. Second, we 
investigate V separately for large and small fe Q ff. Finally, 
we compare our proposed interpolation formula Eq. ([7]) 
to the full range of fe Q ff values. 


A. Dynamics of the model and extraction of the 
long time diffusion constant 

Figures [2] and [3] show the MSD of the tracer particle 
as a function of time, for different unbinding rates k 0 g. 
Symbols represent simulation results. From such plots we 
extract the long time diffusion constant V (see Appendix 
0 by fitting a straight line for large times starting from 
t = r 0 ff (short vertical dashed lines). The results for D 
is shown in Figs. 4][6j but first we discuss some of the 
features of Figs. [2 and [3} 

Figure [2] shows the MSD when crowder particles do not 
diffuse but only unbind and rebind. They rebind either 
always to the same site (upper panel), or to a randomly 
chosen site (lower panel). The short time behaviour in 
both plots is independent of fc 0 ff and is well represented 
by 2Dt (upper dashed dark blue line). The long time 
behaviour is, however, strongly dependent on k 0 g , which 
is evident from the broad scattering of curves. The MSD 
is still linear in time but the diffusion constant (propor¬ 
tional to the extrapolated intersection with the vertical 
axis) depends strongly on k 0 g. The linear regime sets 
in when t ss r a g as is seen from the shorter vertical 
dashed lines. If we increase the particle concentration, 
the shape of the curves remains the same but the scat¬ 
tering of curves increases, since V oc l/p 2 (small k 0 g) 
and V oc 1 — ap (large k 0 g). 

In Fig. [3} crowder particles diffuse and rebind to a 
randomly chose site. As we lower k Q g the separation be¬ 
tween r 0 ff and the collision time t co ii increases, which 
means that the single-file regime ((x 2 (t)} ~ -y/f) becomes 
wider. This is simply because crowder particles have not 
yet started to unbind from the lattice and therefore dif¬ 
fuse collectively in a single-file. We also see that the MSD 
curves for long times is less scattered than before, indi¬ 
cating that V is less sensitive to k 0 g- This agrees with 
our theoretical prediction where V oc y/k Q g compared to 
V oc fc 0 ff when crowder particles stand still. 


B. Small fc 0 ff behaviour 

In Fig. [I] we show how T> depends on small fe 0 ff where 
each panel depicts: (top) immobile crowder particles and 
rebinding to the same lattice site, (middle), immobile 
crowder particles and rebinding to a random lattice site, 
and (bottom) diffusing crowder particles and rebinding 
to a random site. Symbols represent simulation results 
and dashed lines the small k 0 g expressions § ©. Each 
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FIG. 2: Mean squared displacement (x 2 (t)) of the tracer par¬ 
ticle as a function from time for different unbinding rates feoff 
when crowder particles do not move. The crowder particles 
rebind in two ways: to the same site (upper panel), or to a 
randomly chosen site (lower panel). For shorthand we put 
fe 0 ff = feoff/fe c and Toff = Tofffco. Simulation details: lat¬ 
tice constant: a = 1, tracer particle diffusion rate: fc_o = 1 
(ko = 0 for crowder particle), filling fraction: ap = 0.1, num¬ 
ber of lattice sites: 501 (L = 501a), number of simulation 
runs: 9600. 


case is plotted for two concentrations, ap = 0.1 and ap = 
0.2. In order to better compare the three panels with the 
figures below we scaled the vertical axis with the large 
feoff limit D( 1 — ap). In Fig. [ 5 ] we show explicitly how 
this limit is approached. 

The two upper panels in Fig. [4j where the crowder 
particles are immobile are very similar to each other. If 
both cases would be depicted in the same graph, the data 
points would practically sit on top of each other. For 
clarity, we therefore separated the data into two figures. 
We see that the small k Q g behaviour agrees very well with 
the theoretical results, Eqs. (|3| and Q. 
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FIG. 3: Mean squared displacement ( x 2 (t )} of the tracer par¬ 
ticle as a function from time for different unbinding rates k 0 g 
when crowder particles diffuse. Rebinding occurs to a ran¬ 
domly chosen site. For shorthand we put fc 0 ff = fcoff/fco, 
Toff = Tcsko, and r co n = Tcoiifco. Simulation details: lat¬ 
tice constant: a = 1, diffusion rate: kn = 1 (for all particles 
including the tracer), filling fraction: ap = 0.1, number of lat¬ 
tice sites: 501 (L = 501a), number of simulation runs: 9600. 


The lower panel depicts when crowder particles diffuse 
on the lattice. Their movements lead to an overall in¬ 
crease of V for the tracer particle since they no longer 
act as static road blocks. This also changes the scaling 
with k a g from linear in the two upper panels, to Vk 0 g. 
The density dependence is also weaker (1/p compared to 
1 /P 2 )- 


C. Large k a s behaviour 


When k Q ff is much larger than the diffusion rate fco, we 
expect the mean held result Eq. © to hold. We also ex¬ 
pect that corrections to this result becomes increasingly 
prominent as k 0 g is lowered. Both are confirmed by sim¬ 
ulations in Fig. [5j where we see that V/(l — ap) ss 1 for 
k a fi/kD > 1, and D/(l — ap) < 1 for k 0 g/kD < 1. These 
results validate the mean Held argument leading up to 
Eq. § for our quasi one dimensional system. The figure 
only shows the case where the crowder particles rebind 
to the same location, since the behaviour at large k 0 g is 
close to identical for all rebinding modes. 


D. Interpolation formula 

In Sec [HI] we proposed Eq. ([7]) that ties together the 
small and large k 0 g regimes. The comparison to the full 
range of k Q g is shown in Fig. [6] as solid lines (symbols are 
simulation results). Just as in Fig. [4j each panel shows: 
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FIG. 4: Long time diffusion constant T) as a function the 
unbinding rate k D g, when fc D ff is small. Symbols represent 
simulations for the filling fractions ap = 0.1 and ap = 0.2. 
The dashed lines show our predictions Eqs. © <§. Each 
panel depicts: (top) immobile crowder particles with rebind¬ 
ing to the same location (middle) immobile crowder particles 
with rebinding to a random location, (bottom) diffusing crow¬ 
der particles with rebinding to a random location. The data 
points are extracted from linear fits of Figs. [2j|3] (see Ap¬ 
pendix [b]). Jt 2 -values from those fits are larger than 0.98. 
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FIG. 5: Long time diffusion constant T> as a function the 
unbinding rate k a jj, when k 0 a is large. Simulation details are 
the same as in Figs. [2] and [4] 


(top) immobile crowder particles and rebinding to the 
same lattice site, (middle), immobile crowder particles 
and rebinding to a random lattice site, and (bottom) dif¬ 
fusing crowder particles and rebinding to a random site. 
Overall, Eq. |7]) is a good approximation for the whole 
range of fc 0 ff. The deviations are largest in the transi¬ 
tion region, roughly 10~ 3 < k 0 s/kn < 10 _1 , where the 
maximum relative error for all curves is 79% (top panel, 
ap = 0.2). The relative error in the small and large k 0 s 
tails is less than 7%. 


V. SUMMARY AND CONCLUDING REMARKS 

We studied the long-time diffusion constant V of a 
tracer particle in a one dimensional crowded many- 
particle system. We found that V depends strongly on 
the unbinding rate k 0 g of the surrounding crowder parti¬ 
cles and density p. For small k 0 s we made a simple theo¬ 
retical model where we deduced that V ~ k 0 s / p 2 (to first 
order in 1/p 2 ) when crowder particles are immobile and 
only unbind/rebind to the lattice. The prefactor depends 
on how they rebind, either to the same or to a random 
site. When they also diffuse we obtain V ~ yjDk 0 g /p 2 
(to first order in 1/p), a different fc 0 ff-scaling than before; 
D is the free particle diffusion constant. This means that 

V is less sensitive to k 0 s and p when crowder particles 
are diffusing compared to standing still. For large k 0 s, 
we found that all cases agreed with the mean field result 

V ~ D(l — pa), independent of fc 0 ff- Our new expressions 
showed overall good agreement with simulations. 

It is interesting to see which k 0 g / ku regime we expect 
to find in the living cell. As mentioned in the introduc¬ 
tion, residence times of DNA binding proteins vary from 
fractions of a second to up to an hour (unspecific binding 
is even shorter, >5 ms [ 35 ]). To get an order of magni¬ 
tude estimate of k 0 s/k]j, let us assume that fc 0 ff ~ 0.1 
s^ 1 which lies between the in vivo values for the LexA 
and Gal4 transcription factors. One dimensional diffu¬ 
sion constants also have a big variation. They are in the 
range Did ~ 10 5 — 10 7 (bp) 2 s _1 (« 0.01 — 0.1 pm 2 s -1 ) 





Unbinding rate k 0 s/kD 

FIG. 6: Long time diffusion constant T> as a function of the 
unbinding rate k 0 s■ Symbols represent simulations for two 
different filling fractions, ap = 0.1 and ap = 0.2. The solid 
lines shows the interpolation formula Eq. (J7|. Each panel de¬ 
picts: (top) immobile crowder particles with rebinding to the 
same location (middle) immobile crowder particles with re¬ 
binding to a random location, (bottom) diffusing crowder par¬ 
ticles with rebinding to a random location. The data points 
are a compilation from Figs. [4] and [5] 
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[25] , which gives kn = Z?io/(bp) 2 ~ 10 5 — 10 7 s _1 . This 
means that k 0 s/ko ~ 10~ 6 — 10” 8 , which clearly indi¬ 
cates that k 0 s <C kjj. 

The model we studied is inspired by protein diffusion 
on DNA. Our results are simple formulas for the diffu¬ 
sion constant of a tracer particle taking crowding and 
binding/unbinding dynamics into account. Although a 
protein is more complex than a hard-core particle, we 
hope that the simplicity of our results will find its useful¬ 
ness in a range of settings, in particular, single-molecule 
tracking experiments. 
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Appendix A: Simple model for T> when k a s is small 



FIG. 7: Distribution of distances between nearest neighbours 
i p(z) in semi-log scale. Symbols depict simulation results, and 
dashed curves are theoretical results. The filling fraction is 
ap = 0.1. The remaining simulation details are the same as 
in Figs. [2] and [3] 


In this appendix we outline the derivation for the long 
time diffusion constant T> in the small k a g limit that led 
to Eqs. © ©. The main idea is to calculate the typical 
length scale Iq that the tracer travels before being hin¬ 
dered by a crowder particle. In terms of Z 0 , the long time 
diffusion constant is 

V=^, (Al) 

where r is the typical waiting time until a successful 
jumping event. Since the diffusion rate k d is fast, the 
rate limiting step for the tracer particle to move is when 
a flanking crowder particle unbind from the lattice. We 
can therefore envision the tracer particle diffusion as a 
single particle diffusion process on a coarse grained lat¬ 
tice, with lattice constant Iq and jump rate 1/r. 

First we address r, the average time until a successful 
jumping event. Imagine that the tracer particle is flanked 
by two crowder particles, and the time for any of them 
to unbind is 1/2 k Q g. Now, say the that the tracer’s right 
neighbour unbinds but the tracer anyway tries to jump 
left. This jump is forbidden, and so is in fact half of 
all tries the tracer makes. This implies that r is 1 /k a g 
rather than 1/2 k 0 g. Moreover, we also consider rebinding 
of crowder particles, so even if the tracer move in the 
direction of the unbound neighbour it may anyway be 
blocked by another crowder particle. We must therefore 
correct the jumprate with the probability that the site is 
vacant, that is 1 — ap (in equilibrium). In summary, we 
estimate r as 

T = rwr—i ( A2 ) 

fcoff(l - ap) 

Second we turn our attention to the coarse-grained lat¬ 
tice distance Iq- In short, we choose Iq the standard devi¬ 
ation of the distribution of nearest-neighbour distances. 


Here is how we formally arrive at this result. Since it can 
happen that nearest and next-to-nearest neighbours are 
unbound simultaneously, the length that the tracer par¬ 
ticle can move, z, can vary. We choose the probability 
distribution of z to be the probability that there is a sep¬ 
aration 2 between two nearest neighbour particles. The 
distribution of z, ip(z), is known (see Fig. [7|, but differs 
depending on the type of rebinding. If z only can change 
in discrete steps of A, that is A, 2A, 3A..., we can define 
a jump length distribution g(l) for the tracer particle in 
the coarse grained lattice as 

OO 

5 ( 0 = 5(1 + nA)ip(nA), (A3) 


where 6(x) is the Dirac delta function. Now we choose 
lo as the standard deviation of g(l), and from Eq. (A3) 
one can show that 


lo = A \J~(ri% (A4) 

where ( n 2 ) v = n 2 ip(nA). In the subsections below, 

we calculate Iq explicitly for the special cases: (1) immo¬ 
bile crowder particles placed equidistantly, (2) immobile 
crowder particles placed at a random distance apart from 
each other, and (3) diffusing crowder particles. 


1. Case 1: Immobile crowder particles placed 
equidistantly 

In the simulations, the crowder particles sit equidis¬ 
tantly, and unbind and rebind with rates k 0 g and k on , 
respectively. In order to make sure that the average den¬ 
sity p is constant over time, we choose k on = k Q g, and 






work with 2 m particles in the whole system (lattice + 
surrounding bulk). This means that m particles will on 
average be on the lattice and the density is p = m/(aN), 
where N is the number of lattice sites, and a the lattice 
constant. The smallest separation between two crowder 
particles in this setup becomes A = l/(2p), and increase 
discrete steps of A: 


1 13 1 

2 ~p ~p 2 ~p 


(AS) 


Each one of these lengths has a different probability, and 
the distribution of nearest neighbour distances is 





\ (I) 1 "’ 1 ) n = ± 1,±2,... 
0, n = 0, 


(A6) 


which agrees well with simulations (Fig. [7|. Using that 
(n 2 ) v = 6 in Eq. (A4| gives 


l 




(A7) 



FIG. 8: MSD of the tracer particle as a function of time 
for two different unbinding rates. The dashed black lines are 
linear fits of the MSD curves for times t > t 0 s . The simulation 
details are the same as in Fig. [2] 


2. Case 2: Immobile crowder particles with 
rebinding to random locations 

In this case the crowder particles leave the lattice and 
return to a random vacant lattice site. This means that 
the smallest separation is the lattice distance of the orig¬ 
inal lattice, A = a, and distances are in steps of a: 

a , 2a, 3a, ... (A8) 

The inter-particle distance distribution in this case is 


length scale for the coarse grained lattice and a time scale 
associated with a jumping event. 

For this case we cannot use p{z) to estimate Iq since 
<p(z) is the same as when crowder particles are immobile 
(see Fig. [7j O and A), and gives the wrong result for V. 
The reason is that inter-particle distances fluctuate at the 
same rate as the tracer is diffusing, and those fluctuations 
increase T>. In fact, even if k a g = 0, the tracer particles 
still can move across the system, although slowly. We 
estimate Iq as the distance the tracer particle explores in 
a time r, that is 


ip(na) 


^U-ap) 1 " 1 , n = ±l,±2,... 
0, n = 0 


(A9) 


l 0 = \/(x 2 {t = T)) 


/ 4 D 1 -ap \ 1/4 
\irk 0 g p 2 ) 


(All) 


which is corroborated by simulations in Fig. [7J In the 
continuum limit (small a), the distribution becomes ex¬ 
ponential ip(na) ~ e~' n ' ap . Using that (n 2 ) v = (2 — 
ap)/(2a 2 p 2 ), we obtain 



Interestingly, Iq depends on k 0 s and not only p as in the 
previous cases. This changes the scaling of k Q g in V from 
linear in Cases 1 and 2 to y/k a g for this case. This can 
also be understood from the following simple argument. 
The curve for {x 2 (t)) is continuous for all times, and at 
some time the dynamics changes behaviour from single¬ 
file (~ \ft) to regular diffusion (~ t). This occurs around 
t ~ t, which implies 


3. Case 3: Diffusing crowder particles with 
rebinding to random locations 


2 Vt 


t—T 



(A12) 


Here all particles diffuse which drastically changes the 
situation. The main difference is that the tracer does not 
get stuck between two flanking road blocks since they also 
move. However, we know from simulations that the MSD 
for the tracer is in the long time limit proportional to T>t 
(Fig. [3]), which is a direct manifestation that the no¬ 
passing condition is violated (otherwise we would have 
had MSD ~ y/t). Altogether, this implies that there is 


This yields V oc v^off ■ 

Appendix B: Extraction of the long time diffusion 
constant 

The way we determine T> from our MSD simulations, 
is illustrated in Fig. [8] First, T a g is the approximate 
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time at which the MSD becomes linear (shown as vertical 
dashed-dotted lines). Second, we make a linear regression 
of the MSD curve starting from that point, and obtain 
the slope which equals 2V. The resulting fits are shown 
as dashed lines. 


Appendix C: Numerical implementation 

The model (Fig. |T|) is implemented using the Gillespie 
algorithm [40]. The majority of the details of the imple¬ 
mentation has been explained elsewhere [33], but below 
we point out some key differences. 

We keep track of the unbound crowder particles in the 
bulk in order to have the option to rebind them at the 
location they detached from. In practice we use two 


lattices, one of which represents the bulk. The filling 
fraction is maintained at the level we want by setting 
k on = fc 0 ff, and then let the systems equilibrate such that 
half the number of crowder particles sit in the bulk and 
the other half on the lattice we are interested in. This 
representation helpful to investigate all sorts of binding 
modes especially rebinding to the same location. In Ref. 
[34] the bulk served as an infinite particle reservoir and 
the concentration on the lattice was tuned via detailed 
balance (rebinding always occurred to a randomly chosen 
site). Here on the other hand, the bulk has a finite size 
and cannot be seen as a strict particle reservoir. However, 
since we use about 500 particles, fluctuations around the 
filling fraction ap are so small that we rarely (if ever) de¬ 
plete the bulk. This means that we have approximately 
a grand canonical ensemble. 
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